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Abstract 

Usually reason of irreversibility in open quantum-mechanical system is interaction with a thermal 
bath, consisting form infinite number of degrees of freedom. Irreversibility in the system appears 
due to the averaging over all possible realizations of the environment states. But, in case of open 
quantum-mechanical system with few degrees of freedom situation is much more complicated. 
Should one still expect irreversibility, if external perturbation is just an adiabatic force without 
any random features? Problem is not clear yet. This is main question we address in this review 
paper. We prove that key point in the formation of irreversibility in chaotic quantum-mechanical 
systems with few degrees of freedom, is the complicated structure of energy spectrum. We shall 
consider quantum mechanical-system with parametrically dependent energy spectrum. In particu- 
lar, we study energy spectrum of the Mathieu-Schrodinger equation. Structure of the spectrum is 
quite non-trivial, consists from the domains of non-degenerated and degenerated stats, separated 
from each other by branch points. Due to the modulation of the parameter, system will perform 
transitions from one domain to other one. For determination of eigenstates for each domain and 
transition probabilities between them, we utilize methods of abstract algebra. We shall show that 
peculiarity of parametrical dependence of energy terms, leads to the formation of mixed state and 
to the irreversibility, even for small number of levels involved into the process. This last statement 
is important. Meaning is that, we are going to investigate quantum chaos in essentially quantum 
domain. 

In the second part of the paper, we will introduce concept of random quantum phase approxi- 
mation. Then along with the methods of random matrix theory, we will use this assumption, for 
derivation of muster equation in the formal and mathematically strict way. 

Content of the paper is based on our previous studies. However, in this review paper is in- 
cluded also some original material. This part mainly concerns to the discussion about possible 
experimental realization of theoretical concepts in the field of organic chemistry. 

keywords : Statistical Physics, Quantum Chaos, Nonlinear Resonance, Open Quantum- 
Mechanical Systems. 

PACS numbers: 73.23.-b,78.67.-n,72.15.Lh,42.65.Rc 
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Introduction 

The traditional notion of an area, where the laws of statistical physics are effective, 
consists of the assumption that the number of interacting particles is sufficiently large. 
However, a lot of examples of nonlinear systems with a small number of degrees of freedom, 
where chaotic motions occur, had become known by the end of last century 
new stage in the development of notions about chaos and its origin appeared in the last two 
decades of the last century. It turned out that classical Hamiltonian near the separatrix of 
topologically different trajectories may experience a special kind of instability. Because of 
this instability various dynamic characteristics of the system randomly change with time. 
Such a property of the system that performs random motion is called dynamic stochasticity. 
Dynamic stochasticity is an internal property of the system and is not associated with the 
action of some a priory random forces. 

Dynamic stochasticity appears due to the extreme sensitivity of nonlinear system with 
respect to the slightly change of initial conditions or systems parameters. On the other hand, 
even being chaotic, dynamic is still reversible. Irreversibility occurs only after averaging 
of the dynamic over small dispersion of initial data. Note that averaging is not formal 
mathematical procedure. It is essential from the physical point of view. Since dynamical 
description loses its sense due to local instability of phase trajectories. However not the 
existence of initial error is important but, what kind of consequences it has. In case of linear 
system, this influence is negligible. So one can always assume that initial data for linear 
system is defined with the absolute accuracy. But in case of nonlinear systems, even small 
unavoidable error should be taken into account. This leads to the necessity of using concepts 
of statistical physics. As a result, analytical description becomes much more complicated. 
All above mentioned was concerned to the classical case. What is really happening in 
quantum case? Should one still expect non-reversibility in quantum case? The question is 
that as opposite to the classical case, quantum equation of motion is linear. Of course, things 
are more or less clear in case of open quantum systems interacting with the thermostat. If 
so, then, irreversibility appears owing to the averaging of systems dynamics over all possible 
realizations of environment states. Due to this, one can use standard formalism, and from 
the Liouville-von Neumann equation, deduce irreversible in time muster equation, for the 



reduced density matrix 



J]. But how does the irreversibility occur in quantum systems with 



few degrees of freedom This is main question we address in this review paper. 



3 



Usually quantum irreversibility is quantified by the fidelity. Introduced by Peres , this 
concept works pretty well and is especially convenient for the study of problems like quantum 
chaotic billiards [6|. On the other hand, disadvantage is that obtaining of analytical results 
not always is possible and large computational resources usually are needed. 

In this paper, we offer alternative analytical method for the study of problem of quantum 
irreversibility. Our concept is based on the features of the energy spectrum of chaotic 
quantum systems. Peculiarity of energy spectrum of chaotic quantum systems is well-known 
long ago [7|. Random Matrix Theory (RMT) presumes eigenvalues of chaotic systems to 
be randomly distributed. Number of levels also should be quite large. In order to deduce 
muster equation for time dependent chaotic system, we will utilize methods of RMT in the 
second part of this paper. However, what we want to discuses in the first part is completely 
different. The key point is that, study of quantum chaos as a rule is focused on the semi- 
classical domain. We mean not only Gutzwiller's semi-classical path integration method j^, 
but also RMT. Since the RMT, in somehow implies semi-classical limit. At least implicitly, 
due to the large number of levels included into process. 

In the first part of our paper, we shall consider chaotic quantum-mechanical system with 
few levels included into process. In spite of this, feature of the energy levels leads to the 
irreversibility. 

Paper is organized as follows: 

In the first part we shall consider quantum mechanical-system with parametrically de- 
pendent energy spectrum. This parametrical dependence is quite non-trivial, contains do- 
mains of non-degenerated and degenerated stats separated from each other by branch points. 
Namely energy spectrum of our system is given in terms of Mathieu characteristics. Due to 
the modulation in time of the parameter, system will perform transitions from one domain 
to other one. For determination of eigenstates for each domain and transition probabilities 
between them, we utilize methods of abstract algebra. We shall show that peculiarity of 
parametrical dependence of energy terms, leads to the formation of mixed state and to the 
irreversibility, even for small number of levels involved into the process. So this part may be 
considered as an attempt to study quantum chaos in the essentially quantum domain. This 
study is based on our previous papers j^, [lo|, ll|, Q • In addition, in the present paper, we 
will discuss in details possible experimental realizations and applications of the theoretical 
concepts in the organic chemistry and polyatomic organic molecules. 
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In the second part of the paper, we will introduce concept of random quantum phase 
approximation. Then we will use this assumption, for derivation of muster equation in the 
most formal and mathematically strict way. 

I. QUANTUM PENDULUM 
A. Universal Hamiltonian 

Let us present the atom as a nonlinear oscillator under the action of the variable 
monochromatic field. Then the Hamiltonian of the system atom + field is of the form: 

H{x,p,t) = Ho{x,p) + Hnl{x) + eV{x,t), (1) 

where 

2 

i/o = l/2(^ + cUoW), Hnl= Px^ + ix^ + (2) 

V(x,t) = Voxcosnt, eVo = — /o, e < 1. (3) 

m 

Here x and p are coordinate and the impulse of the particle (electron), uq is the frequency of 
oscillations, /5 and 7 are coefficients of nonlinearity, m and e are the mass and charge 
of the particle, /o is the amplitude of the variable field. Having made passage to the 
variables of action-angle 1,6 with the help of transformation x = {21 /mtooY^'^ cos9,p — 
— (2/mu;o)^^^ sin 6*, supposing resonance condition f2 t^o is realized and averaging with 
respect to the fast phase 6, one can obtain: 

H = Ho{I)+eV{I) cosv^, 

where 

Ho{I) = ivol + H^L, H^L{I) = ^{^ri (4) 
ip = e-ujt, eV{I) = ^/l/muoVo. (5) 

The role of nonlinear frequency plays (dHo/dl) = uq + ujnl{I) where unl{I) = 
(37r/2)(/7/mci;Q). Now suppose that nonlinear resonance condition luq + loj^tl^Iq) = f2 is 
fulfilled for action / = Iq. It is easy to show that for a small deviation of action from the 
resonance value A/ = I — Iq (A/ ^ Iq) after the power series expansion, if the condition of 
moderate nonlinearity is just /i <^ 1/e, where 

{dUNL{I)/dI)\l 
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for Hamiltonian we obtain: 

H='^{Alf + eV{Io)cos^. (6) 

As usual (6) is called universal Hamiltonian. Let us notice that Hamiltonian H is the 
Hamiltonian of pendulum, where I jio^i^ plays the role of mass, A/ plays the role of the pulse 
and V{Iq) plays the role of potential energy. If in (6) A/ is substituted by the appropriate 
operator A/ —ifid/dip, one can obtain the universal Hamiltonian in the quantum form 

With the help of (7) it is possible to explore quantum properties of motion for the nonlinear 
resonance. 

Having written the stationary Schrodinger equations 



H^n = En^jjn (8) 



for the Hamiltonian (7), we get 



+ (E„-\/(/,(^))^/;„ = 0, (9) 
V{l,if) = 2/cos2(/j 



where the dimensionless quantities are introduced 



^..-Ife.'-^ (10) 



and the replacement (/) — > 2(y9 is done. 

The interaction has the following properties of the symmetry: \.V{ip) = 
V{-ip), 2.V{ip) = V{7r + if), 3.V{^ + = - ip). G.M.Zaslavsky and G.P.Berman 
[8] were the first who considered the equation of the Mathieu-Schrodinger for the quantum 
description of the nonlinear resonance in the approximation of moderate nonlinearity. They 
studied the case of quasi-classical approximation A/ ^ h for both variables /q and A/. In 
this review we investigate the equation of the Mathieu - Schrodinger in essentially quantum 
area. 
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FIG. 1: The dependence of the energy of interaction V from the phase. 
B. Periodical solution of Mathieu - Schrodinger equations 

We content ourselves with only even and odd solutions with respect to ip of equation of 
the Mathieu - Schrodinger (9). Those solutions have n zeroes in the interval < (f < n. 
Eigenfunctions ipn can be recorded with the help of the Mathieu functions [9,10]: even 
ce„(v9, /) and odd se„(v9, /). Appropriate eigenvalues are usually designated by a„(/) and b{l). 
For simplicity below sometimes we omit the argument / and write ce„ se„ (</)), a„, bn- 
Mathieu functions are eigenfunctions of the problem of Sturm-Liouville for the equation (9) 
for the boundary conditions 



dip dip 



0, for se„(v2) 
= 0, for ce„ ((/)). 



(11) 



From the general theory of the Sturm-Liouville it follows, that for arbitrary n = 1,2,... 
there exists eigenfunction scnif) and for each n = 0,1,2,... determined eigenfunction 
ce„ ((/:>). The definition of the Mathieu function must be supplemented with the choice of the 
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arbitrary constant so that the conditions were fulfilled: 

cen(0,/)>0, -/ cel{Lp,l)d(f 



IT 







-(0,/)>0, - / sel{ifj)dif = l. (12) 

Jo 



dsen , ' '"^^ 

dip 

If ip{ip) = G{ip) means either ce„(v9) or se„,(v9), then G{ip) and G'(7r — ip) satisfy the same 
equation (9) and the same boundary conditions (11). Therefore these functions differ from 
each other only by the constants. Hence, G{p) is even or odd function with respect to 
11/2 — p. Taking this into account, two functions (11) break up into four Mathieu functions: 

t/;(0) = i){-K/2) = 0, G{ip) = Se2m+2{^), phttSC 71, 

^(0) = ^(7r/2) = 0, G{p) = se2ra+i{'f), phase 2%, 
dp 

^(0) = V^(7r/2) = 0, G{p) = ce2m+i{p), phase 27r, 
dp 

3^(0) = 3^(^/2) = 0, G{p) = ce2mip), phase vr. (13) 
dp dp 

For arbitrary m = 0, 1, 2, . . . there is one eigenfunction for each of four boundary condi- 
tions and m equal to the number of zeroes in the interval < p < n /2. The functions (13) 
represent a complete system of eigenfunctions of the equation (9). 



C. Symmetries of the equations of Mathieu-Schrodinger 

The properties of symmetry of the Mathieu-Schrodinger equation can be presented in 



Table 1 



15|- 



As is known, group theory makes it possible to find important consequences, following 
from the symmetry of the object under study. Below, with the aid of group theory we will 
establish the presence (or absence) of degeneracy in the eigen spectrum of the Mathieu- 
Schrodinger equation and a form of corresponding wave functions By immediate check it is 
easy to convince, that four elements of transformation 

G{p —p) = a, G{p ^ TT — p) = b, 

G{p ^ IT + p) = c, G{p —></)) = e 
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Gi<p) 


G{-9) 


G{7T — (p) 


G{;r + <p) 




ce,„{9) 


CBl^ {9) 


{9) 




ce:,,.! {9) 


-csj^^i {9) 


-ce.^^i {9) 


SSj^.l {<p) 


-i-e^M+i {9) 


sei^.i {9) 


-ss._^.A9} 


sei^.i {9) 


-sej^^j {9} 


-SS2^A9) 


se2^.2{9) 



TABLE I: Relations of a symmetry for the Mathieu function. 

are forming a group. For this it is enough to test the reahzation of the following relations: 

2 1,2 2 
a = = c = e, 



ab = c, ac = b, be = a. (14) 

The group contains three elements a, b, c of the second order and unity element e. The 
group G is isomorphic to the well-known quadruple group of the Klein ly, [iTj. This group 
is known in the group theory by the applications to the quantum mechanics (designated as 
V^). All the elements of the group commute. This assertion can be easily checked taking into 
account group operations (14). So, the symmetry group of the Mathieu function G is the 
Abelian group and consequently has only one-dimensional indecomposable representations. 

The presence of only one-dimensional representations of the symmetry group, describ- 
ing the considered problem, hints on the absence of degeneration in the energy spectrum. 
So, we conclude, that the eigenvalues of the equation of the Mathieu - Schrodinger (9) are 
non-degenerated, and the eigenf unctions are the Mathieu functions (13). However we shall 
remind, that both the energy terms a„, and the Mathieu functions depend on the pa- 
rameter /. At the variation of / in the system can appear symmetry higher, than assigned 




FIG. 2: Four vectors that are transformed by the elements of quadruple group V. 

in Table 1, that might lead to the degeneration of levels. 

In order to make more obvious the isomorphism of the symmetry group G{lp) of Mathieu- 
Schrodinger equation with quadruple Klein group V, let us consider the plane of rotation of 
phase ip. 

The orientation of vector v4^„^, presented in Fig. 2, is obtained from by means of mirror 
reflection relative to the plane passing through the axis Oy perpendicular to the figure plane 
ay. The orientation of vector v4^_|_<^ is obtained by means of rotation by angle vr about an 
axis, passing perpendicular to the figure through the origin of the coordinate system C2 (axis 
of rotation of the second order). The symmetry elements a^, ay and together with the unit 
element form quadruple group V : e, a^, ay, €2- Now it is possible to bring the elements 
of two group to one-to-one correspondence: G{—ip) ax, G{tx — ip) ^ ay, G{tx + </?)—>■ 
G2, e — e that proves the isomorphism of above mentioned groups. 

Each of three elements a, b, c in combination with unit e forms a subgroup 

D G+ : e, b, 

G D : e, c, (15) 
D Go : e, a. 

Each of subgroups Gq, G+, G_ is the invariant subgroup in-group G. The presence of the 
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FIG. 3: Fragment of energy spectrum of quantum pendulum, composed from theoretical-group 
consideration. 

three invariant subgroups of the second order indicates the existence of three factor-group 
(4:2=2) of the second order. 

Fo : E(e,a),A(6,c), 
F+:E(e,6),A(a,c), 

F_ : E(e,c),A(a,6). (16) 

Here we introduced notations, common for the group theory: E stands for a unit element 
of factor-group; A stands for an element of factor-group. Group G is homomorphous to its 
factor-groups Fq and F±. As we see from (16) the elements of factor-group are formed as 
a result of unification of certain two elements of group G. This kind of unification of the 
elements of group G indicates pair merging of energy levels ce2m, ce2m+i, 562™+! and se^m^i- 
Therefore, three kinds of double degeneration, corresponding to the three factor-groups Fq 
and -F-J-, appear in the energy spectrum. This happens owing to the presence of parameter 
/. Therefore, for the different values of parameter / energy spectrum may be different at 
least qualitatively. Let us find out how to present combination of such a variety in energy 
spectrum 

Factor-group Fq is responsible for pair unity of levels of the same symmetry relative to 
the center of potential well G{t\ — ip) (Fig.l). Quantum vibration motion appears because 
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FIG. 4: The group of the eigenvalues of a„(/) and h^{C) as functions of I are plotted by numerical 
methods jisl - 

of this unification. In Fig. 3. these levels are arranged to the extreme left and right from the 
ordinate axis. 

Factor-groups Fq and are responsible for pair unity of equations ce„ and se„(?2 = 1,2,...) 
and for formation of clockwise and counter clockwise quantum rotation motion. In Fig. 3 
corresponding degenerated levels are found in both sides close to the ordinate axis. Numer- 



ical calculations 



141 ] (see Fig. 4), of Mathieu characteristics aiX) and hiV) prove that energy 



spectrum of quantum pendulum has a very complicated form and manifests all the main 
features of the spectrum presented in Fig. 3. 

Similarity of the plots presented in Fig. 3 and Fig.4 is obvious. Namely, both of them hold 
symmetry with respect to Y-axis. From both sides of Y-axis same eigenstates are coupled 
(degenerated). Only difference is that Fig. 3 belongs to the theoretical group analysis, while 
Fig.4 is plotted using numerical methods. Both of these methods have their advantage and 
disadvantage. Only way to evaluate exact positions of branch points is to use numerical 
methods. From the other hand, eigenfunctions for each domain should be defined using 
theoretical algebraic methods. For more details see lOj. 
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D. The Physical Problems are Reduced to Quantum Pendulum 



We have already met with one physical problem that can get reduced to the solution of a 
quantum pendulum (9) - this is problem of quantum nonlinear resonance. Now we will get 
to know with other quantum-mechanical problems, that also get reduced to the solution of 
quantum pendulum. 

As is known jisl . one of the forms of internal motion in polyatomic molecules is 
torsion oscillation which for sufficiently large amplitudes transforms to rotational motion. 
In order to describe the corresponding motion in Hamiltonian we assume that ip is the angle 
of torsion of one part of the molecule with respect to the other part and replace the mass 
m by the reduced moment of inertia / = /i/2/(/i + I2), where Ji and I2 are the inertia 
moments of rotation of the parts of the molecule with respect to its symmetry axis. Thus 
we obtain 12| 



f/(^) = ^(l-cosn^) (17) 

where Vq defines the height of potential barrier that separates torsion oscillations from the 
rotation of one part of the molecule with respect to the other part, and n defines the quantity 
of equilibrium orientations of one part of the molecule with respect to the other part. For 
the molecule of ethane H^C — CH^ , dimethylacetylene H-^C — C = C — CH^ and for other 
organic molecules we have n = 3 equilibrium configurations (see Fig.5.). 

The configuration shown in Fig.5 corresponds to an energy maximum and is a 
non-equilibrium configuration (cis-configuration). Other non-equilibrium configurations 
are obtained by rotating by the angles ^ and 2^. Equilibrium configurations (trans- 



configurations) are obtained by rotating o: 
Below we give the numerical values pi 



the angles f , ^ + ^, f + 2^. 



3 ' 3 ' 3 



191 of other parameters of some organic molecules 



having the property of internal rotation. Thus for the molecule of ethane C2HQ we have 
= I2 ^ 5.3-10"'^'^ kg ■m'^, Vo{C2Hq) ^ 2. 1-10^^° J, and for the molecule of dimethylacetylene 

C^He we have h = h ^ 10.6 ■ 10-"^^% ■ m\ VoiC^He) ^ 0.34 ■ lO'^oj. 

The Schrodinger equation corresponding to Hamiltonian (7) has the form 

+ ft^^ - ^^0(1 - cosn^)]^ = 0, (18) 

where Ek is the eigenenergy of the k — th state. Note that Sk = £k{Vo) is the function of 
barrier height Vq. The condition of motion near the separatrix (near a potential maximum) 
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FIG. 5: A schematic drawing of the molecular structure of ethane H-^C — CH-^. The circular arrow 
shows the torsion phase 93, tq is the equilibrium distance between two parts of molecule. 




33 33 



FIG. 6: Potential energy curve of torsion motions in H3C — CH3. Phase ip = corresponds to 
equilibrium trance-configuration. 

is written in the form Ek ~ K. If we introduce the new variable a = then equation (18) 
can be rewritten as 

— + [E- 2/0 cos 2a]tlj{a) = 0, (19) 
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FIG. 7: a) The relative motion of atom Li to hard fragment C = N. b) Coordinates r and (p 
describe relative motion Li to fragment C = N. 

where 

E = ^^ie-Vo/2) (20) 
plays the role of energy in dimensionless units, and the parameter 

lo = 44^0 (21) 

is the half-height of the barrier in dimensionless units and plays the same role as the length 
of the thread does in the classical pendulum problem. 

As our next example, we consider the vibration dynamics of a triatomic floppy molecule: 



the LiNC/LiCN isomerizing system which has been extensively studied 20|, l2l|. This 
molecule presents two stable isomers corresponding to the linear configurations, Li — NC 
and Li — CN, which are separated by a relatively modest energy barrier. The motion in the 
beginning is very floppy, and then the Li atom can easily rotate around the CN fragment 

One can describe relative motion of Li to hard fragment C = by means of two coor- 
dinates: (p - angle of orientation of Li with respect to hard fragments axis and r distance 
between atom Li and mass center of fragment C = N. Angle ip = corresponds to Li — CN 
and ip> = n corresponds to its isomer CNLi. It is easy to make sure, that such isomerization 
process may be described by using the potential (17), when n = 2. We obtain Schrodinger 
equation for such case in the same form, as we had for previous one (19). 
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E. Degenerate states of Mathieu-Schrodinger equation. 



In the theory of Mathieu functions, the graphs of the eigenvalues a„(/) and bn{l) as func 



tions of / are plotted by numerical methods 



151 ] ■ As seen from these graphs, curves a.„(/) 



and &n(0 merge for small /, while curves a„(/) and bn+i{l) merge for large /. It is obvi- 
ous that the merged segments of the Mathieu characteristics correspond to the degenerate 
states whose existence has been mentioned above. In this section, we will define the wave 
functions of degenerate states. Below, the presence of branch points will play an essential 
role in explaining the transition from the pure state to the mixed one during the quantum 
investigation of the dynamics near the classical separatrix. In what follows, we will use the 
plane with coordinates /, E. In the classical consideration, the motion of a mathematical 
pendulum in a neighborhood of the separatrix occurs when the initial kinetic energy of the 
pendulum is close to the maximal potential one. It is obvious that, on the plane {l,E), to 
this condition there corresponds the straight line E = 21. Therefore we can say that, on 
the plane (Z, -E), to non-degenerate states there corresponds a certain domain lying on both 
sides of the line E = 21. It is in this very domain of the change of / that the system is 
characterized by symmetry group G. 

a) Degeneration of states at small I {area from the left of the separatrix line) 

In the limit / — > the equation of the Mathieu - Schrodinger (9) takes the form: 

^ + i^n^n = (22) 

The orthonormal system of solutions of the equation (22) consists of even and odd solutions 

ifjg = COS nip, tjju = sin nip. (23) 

They both correspond to the same energy value E^ = n^. Note that functions (23) corre- 
spond also to the well-known asymptotic (/ 0) forms of the Mathieu functions [9]: 

cen{ip) — ^ cosnip), sen{ip) — > sm.nip. (24) 

This means that at the diminution of / the coming together of the energy terms with the 
identical n takes place and for / = they are merged together. It is necessary to find out 
that this confluence happens at the point / = or at / = Z^'' ^ 0. In this section, below we 
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will be concerned with finding a lower point of the merging of terms Z^''. At first let us find 
out what the eigenfunctions of the degenerated states corresponding to the level En = v? 
look like. Equation (22) is the Schrodinger equation for free rotation in the phase plane (/j. 



The continuous Abelian group of two-dimensional rotations 0^(2) IGj corresponds to this 
motion. 

Since the Abelian group may have only one- dimensional irreducible representations, the 
two-dimensional representation constructed in the base of real-valued functions (23) will be 
reducible. Hence functions (23) cannot be eigenfunctions of a degenerate state. To surmount 
this problem we shall recollect that the eigenfunctions for the degenerate condition can be 
also complex. 

As is known symmetry relative to the change of time sign in the Schrodinger equation, 
accounts for the fact that the complex-conjugate wave functions correspond to one and the 
same energy eigenvalue. Therefore two complex-conjugate representations ipn{'^) and ipniv^) 
should be regarded as a representation of doubled dimension. Usually, for the basis of the 



irreducible representation of the group 0^(2) complex functions are assumed [161] 



^|J^{ip) = e-^"^. (25) 

Therefore, in the degenerate domain, in view of conditions of normalization, following com- 
plex conjugate functions should be considered as eigenfunctions 

MV) = ^ rniv) = ^e-^. (26) 

Let us remark that group 0^{2) is isomorphic to subgroup (24). The element of the 
symmetry c = G{ip ^ vr + (/?) of the subgroup provides recurrence of the phase variation 
after each period and consequently the symmetry G- characterizes the condition of motion 
similar to the classical rotary motion. However, to use only the argument of symmetry is 
not sufficient for finding the coordinates of the branching point Below to find these 

points we use the secular perturbation theory. So, at / = we have doubly degenerate states 
with the wave functions (26). Let us find out, whether the perturbation 

V{l,ip) =lcos2ip, /<1 (27) 

can remove the existing degeneration. It is known that, first order terms of the perturbation 
theory for the energy eigenvalues and the exact functions of zero approximation for double 
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degenerate levels look like [1 



4± = ^[(^11 + V22) ± V(l^ll + V^22)2 + 4-|V] 



(0) 



'12 



(0) 



± 



2\Vu\ 
V21 



(0) 

1 ± ■ 



1^11 - V, 



'22 



121^], 



1/2 



21V, 



12 



It 



^22)2 

Vn 



4. |y.J2 



12 J 



22 



1/2 



(28) 



^J{Vll + V22f+^-\Vi2\\ 

where the index in brackets corresponds to the order of the perturbation theory, matrix 
elements of the perturbation (27) Vik{i,k = 1,2) are calculated by using of functions (26) 
of the degenerate state of the non-perturbed Hamiltonian. Taking into account expressions 
(26) we shall calculate the matrix elements: 



Vu = I I i^nif) ■ i^nif) cos 2ipdip 

'0 



0, V, 



22 



0, 



Vi 



12 



TT 



COS 2(pdip 



if n ^ 1 
^ if n = 1 

After substitution of those matrix elements in the expressions (28) for the eigenvalues and 
exact eigenf unctions we shall obtain: 

Itt 



E. 



(1) 



V'n=i = cos(^, V'„=i = -^sin(^. 



(29) 



Thus, the exact wave functions (29) of the non-degenerate state only for = 1 coincide with 
the Mathieu function in the limit / (24). 

The perturbation V{l,ip) removes degeneration only for the state n=l. Therefore it is 
only for the state n = 1 that the spectrum branching occurs at the point / = 0, which 
agrees with numerical calculations given in the form of diagrams (see Fig. 4.). It can be 
assumed that in the case of diminishing /, the merging of energy terms for states n ^ 1 
takes place at the point at which the states are still defined by the Mathieu functions and 
not by their limiting values (24). Wave functions for degenerate states / 7^ 0, n ^ 1 can 
be composed from the Mathieu functions by using the same arguments as have been used 
above in composing the wave functions for / 0). As a result, we obtain 

V2, 



2m+l 



^n(/,V5) 



-(ce„v9 ± isen^), n = 2m + 1, 
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FIG. 8: The energy levels as function of parameter / on the plane {E, I) on the left area from the 
separatrix line. The points of the branching of curves represent the boundaries between degenerate 
and non-degenerate states. 

^ ^ ^{cen^ ± isen^), n^l, 1^0, n = 2m + 2. (30) 

Let us assume that at /„ = the removal of degeneration for the n — th energy term 
happens. 



b) Degenerate state at major I. The area on the right of the separatrix line. 

With the increasing of / the particle can be trapped in a deep potential well {V = 
/cos 2^9, < ip < n Fig. 1.), and perform oscillatory motion. Properties of wave functions 
of quantum oscillator near to the bottom of the well are well known. This is the alternation 
of even and odd wave functions relative to the center of the potential well 7r/2 and presence 
of zeros in wave functions. With the help of the third column of Table 1 it is possible to 
write symmetry conditions close to ir /2: 

71 TT 

cem{^ + f) = (-l)'"ce^(- - ^) 
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sem{- + v) = (-1) 



m+l 




(31) 



i.e. ce2m{f), se2m+i{f) are even functions and se2m{f), ce2m+i(v^) are odd functions. Func- 
tions ce2m(v^), se2m+i{f), ce2m+i{f) and se2m+i{f) have m real zeros between ip = and 
if = n/2 (not considering zeros on edges). The existing alternation of states (Fig. 8.) in area 
along the line of the separatrix is conditioned by the properties of states at the small /. With 
the help of the expressions (31) it is possible to determine easily, that in the spectrum of the 
states along the line E = 21 two (instead of one) even states alternate with odd states and so 
on. To get the alternation, caused now by properties at major /, two even conditions must 
degenerate in one even and two odd - in one odd. So we come to the conclusion, that two 
levels with wave functions ce2m(v') and se2m+i{f) coming nearer amalgamate in one level, 
and the following two levels ce2m(v^) and se2m+i{'^) also in one level. The levels obtained in 
this way will be doubly degenerated. It can be assumed that with the growth of / the states 
defined by the symmetry group G transform to the states with the symmetry of an invariant 
subgroup (15). This transformation takes place at the merging point of non-degenerate 
terms /.„ = Z^''. Recall that subgroup G+ contains two elements: the unit element e and 
the reflection element with respect to the symmetry center of the well b = G{ip ^ n — ip). 
Complex wave functions of the area of degenerate states, with the symmetry of the invariant 
subgroup G+, can be composed of pairs of functions of merged states in the same manner 
as we have done above for the area of small I for states with the symmetry of G^. 

Not iterating these reasons, we shall write complex wave functions corresponding to the 
degenerated states in the form 



In the base of complex wave functions ^2n and C2n+i ^he indecomposable representation of 
the subgroup G+ (15), is realized. Evenness of the wave functions ^ and ( with respect to the 
transformation b = G{(p tt — (f) of the subgroups G+ characterizes an important property 
of wave functions evenness of the quantum oscillatory process. The results, obtained in this 
section, are plotted in Fig. 9. 

Figures 8 and 9 supplement each other: in the field of intersection with the separatrix the 
curves of the Figs. 8 and 9 are smoothly joined. So, we shall add up outcomes obtained in 



^2n{'^) = ce2n{<^) ± se2n+i{<^) evcu statc 



C2n+M = ce2n+i(v5) ± se2n+2{^) odd stttte 



(32) 
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FIG. 9: Energy levels as a function of the parameter / on the plane {E, I) to the area right from 
the separatrix line. The points of branching of curves represent degeneration points of terms in 
this area. 

this section. The Mathieu-Schrodinger equation has an appointed symmetry. The transfor- 
mations of the symmetry of the Mathieu functions form group G, which is isomorphic to the 
quaternary group of Klein. To this symmetry on a plane {E, I) corresponds the appointed 
area along the line of the separatrix E = 21, containing non-degenerated energy terms. This 
area is restricted double sided by the areas of degenerate states, which are characterized by 
the symmetry properties of the invariant subgroups G- and G+, respectively. The bound- 
aries of these defined by the branching points of energy terms existing both on the 
right and on the left of the separatrix. 

The area of degenerate states is the quantum-mechanical analogs of two forms of motion of 
the classical mathematical pendulum-rotary and oscillatory. Comparing results of quantum 
reviews with classical, we remark that these two conditions of motion at quantum reviewing 
are divided by the area of a finite measure, whereas at the classical reviewing measure the 
separatrix is equal to zero. 
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F. Quantum analog of the stochastic layer. 



In the case of Hamiltonian systems, performing a finite motion, a stochastic layer formed 
in a neighborhood of the separatrix under the action of an arbitrary periodic perturbation is 
a minimal phase space cell that contains the features of stochasticity |l| . In this section we 
shall try to find out what can be considered as the quantum analog of the stochastic layer. 

Let us assume that the pumping amplitude is modulated by the slow variable electro- 
magnetic field. The influence of modulation is possible to take into account by means of 
such replacement in the Mathieu-Schrodinger equation (9), 

I lo + Alcosut, Al<lo. (33) 

Here A/ stands for the amplitude of modulation in dimensionless unit (see (10)), u is the 
frequency of modulation. We suppose, that the slow variation of 1 can embrace some quantity 
of the branching points on the left and on the right of separatrix line (Figs. 8,9) 

A/> \ll-r\, n = l,2,...N. (34) 

As a result of replacement (33) in the Hamiltonian (7), we get 

H = Ho + H'{t), (35) 

H'{t) = Al cos 2ip cos ut, (36) 

where Hq is the universal Hamiltonian (7) and H'{t) is the perturbation appearing as a 
consequence of pumping modulation. 

It is easy to see, that the matrix elements of perturbation (36) H'{t) for non-degenerate 
states equal zero. Really, having applied expansion formulas of the Mathieu functions in the 
Fourier series [l^ it is possible to show 

l'2lT 

< cen\H'{t)\sen >^ / ccni^p) cos2Lpse2n{^)d(p = (37) 

simultaneously for the even and odd n. The expressions of the selection rules (37) will 
be fulfilled for values / from the area /" < / / Transitions between levels cannot be 
conditioned by time-dependent perturbation (36). It is expedient to include perturbation in 
the unperturbed part of the Hamiltonian. The Hamiltonian, obtained in such way, is slowly 
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depending on the parameter /. So, instead of (35) and (36) for the non-degenerated area G 
we get the Hamiltonian in the form 

H = —^ + l{t)cos2^, (38) 

l{t) = lo + Al cos ut. (39) 

There arises the situation in which the system slowly "creeps" along the Mathieu charac- 
teristics and, in doing so, encloses the branching points on the left /" or on the right 
a) Irreversible "creeping" of energy term populations due to the influence of a 
measuring arrangement. 

According to the general rules of quantum mechanics, probabilities that the system will 
pass to the eigenstate of another area are defined by the coefficients of expansion of the wave 
function of one area into the eigenfunctions of another area. Let us assume that, initially, 
the system was in one of the eigenstates from the non-degenerate area G, for example, in 
the state ce2n- After a quarter of the modulation period T/4 (where T = 2p/n), having 
passed through the point the system finds itself in the degenerated area G^. In this case 
the system will pass to degenerate states with probabilities, 

■ vr 



P(ce2„ ^ iptJ = 1-1 ce2ni^)ipi*{(f)dip\' 







1 f 

= — I / ce2n{f){ce2n{f) ±ise2n{f)Tdip\^ = 1/2. (40) 
27r Jo 

For deriving (40) we used the condition of normalization (12) and orthogonality [9] 

/■2n 

/ cefc(v9)sei+i((/?)rf</? = 0, /,/c = 0,1,2, ...). (41) 
Jo 

The passage (40) is based on the assumption of having a deep physical sense. As is 
generally known, in quantum mechanics symmetry with respect to both directions of time 
is expressed in the invariance of the Schrodinger equation with respect to the variation of 
the sign of time t and simultaneous replacement tp by way of tp*. However, it is necessary 
to remember that this symmetry concerns only the equations, but not the concept of a 



measurement playing a fundamental role in the quantum mechanics |l7l.l22[|. "Measurement" 
is understood as the process of interaction of the quantum system with the classical object 
usually called "instrument." Under the measuring arrangement, consisting of the analyzer 
and detector, one must not imagine the laboratory's instrument. In our case, the role of the 
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instrument plays in our case the modulating field, which is capable to "drag" the system 
through the branching points. When passing through the branching point from one area 
to another, the state remains unchanged. However, being an eigenstate in one area, it will 
not be an eigenstate in another. At the passage through branching points there occurs a 
spectral expansion of the initial wave function belonging to the region of one symmetry 
over the eigenf unctions belonging to the region of another symmetry. The presence only of 
the analyzer reserves a pure state and the process remains reversible. So, passage through 
the branching point plays role of analyzer. Further we shall assume the presence of the 
detector, defining which of the states ip^ V'n is involved in passage. The transition of 
the system to various states defined by probabilities (40) is fixed by means of the action of 
the detector. The presence of the detector is expressed formally in averaging with respect 
to phase and neglecting the interference term usually appearing in the expression for a 
distribution function. As a result of averaging the partial loss of information about the 
condition of the system takes place and a mixed state is generated. In our problem role 
of detector will play effect self-chaotization, which appeared in degenerate region (see next 
subsection) . 

As is follows from (40), after the quarter period degenerated rotary states ^/^^ and il)2n 
will be occupied with the identical probability. After the half period |T the system again 
appears in the area G going through the branching point /" in the reverse direction. At the 
same time there appear probabilities, of the transition into the states ce„, se„, and both of 
them are distinct from zero 

1 1 f"^"^ 1 
Pi-i^tn ce2„) = 7tI- / (ce2n(^) ± ise2n{^))ce2n{^)d^? = -, (42) 

1 1 f^'^ 1 

-P(V'2n ^ ■se2n) = ttI- / ice2ni'^) ± ise2n{^))se2ni^)d(p\'^ = -. (43) 

2 TT Jo 2 

Here we have used again normalization (12) and orthogonality relations (41). It is easy to 
write transition probability from ce2n into one of the degenerated states and back in the 

ce2n 

P{ce2n ^ Ce2n) = P{ce2n ^ ^ Ce2n) = 

= P{ce2n i^L)Pi^L Ce2n) + P(ce2„ ^ i^2n)Pi^2n ^ Ce2„). (44) 

Here the first summand corresponds to the passage through the degenerated state ip2n ^-nd 
the second one to the passage through It is easy to see with the help of previous 

24 



computations (40), (42), and (43) that contributions of these passages are identical and 
individually equal to 1/4. Therefore finally we have 

P_(ce2n ^ ce2„) = ^. (45) 

Similarly it may be shown that transition probability from the state ce2n in one of the 
degenerated states t/'^ and back in the area G, in the state se2n by means of going through 
the point Z" is 

P{ce2n ^ Se2n) = P{ce2n ^ '4'2n)P{'^2n ^ Se2n) + 
+ P(ce2„ -> i^2n)P{i^2n Se2n) = \-\ + \-\ = (46) 

Thus, the system being at the initial moment in the eigenstate ce2n, at the end of half-period 

of modulation appears in the mixed state p2n in which the states ce2n and se2n are intermixed 

with identical weight, and corresponding levels are populated with identical probabilities. 

After the expiration of quarter of cycle the system will pass from the area G (the state 

p2n) in the area 6*+, going through the point /+. In passages from the area G+ four states 

take part S,fn = ■^{c^2n ± ^■se2n+i) and Ctn-i = "^(ce2n-i ± ise2n)- So, with taking into 

consideration the above mentioned for the probabilities of transitions we get 

1 f'^'^ 1 
Pip2n ^ ^tn) = -r\ / i^2n{^) + se2„ (</?)) (ce2„ ((p) T ise2„+i (^))c/</?P = -, (47) 
47r Jo 4 

1 /"^'^ 1 

P{p2n Ctn) = / ice2n{^) + Se2ni^)){ce2n^l{^)Tise2ni^))d(p\'^ = (48) 

47r Jo 4 
For deriving the last expressions in addition to the normalization conditions we have used 

the orthogonality conditions 

r2-K /•2tt 

/ cen{(p)cemi'p)d(f = / sen+MsCm+Mdip = 0, m^n. (49) 
Jo Jo 

On the basis of (47) and (48) we conclude, that after the time |T, system will be in the 

area G+ in one of four oscillatory states .^^„i and C2n-i with the identical probability equal 

to 1/4. 

After one cycle T the system gets back in the area G, from which it started transition from 
the level ce2n- Upon returning, four levels ce2n, se2n, ce2n-i and se2n+i will be involved. 
Calculating probabilities of passages from the oscillatory state of the area G+, to these four 
levels we shall obtain 

Pi^tn - Ce2n) = P{^tn ^ ^62^+1 = 1/2, (50) 
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^(^2^1 Se2n) = P(C2n-l ^ Se2n^, = 1/2. (51) 

The probability of passages from the nondegenerated area to the area G+ in one of the 
oscillatory states ^f^, Ctn-i back in the area G will be: 

P+{p2n ^ Se2n) = P{p2n ^ ^2n)P{^2n ^ Se2„+l) + 
+ P(p2n - ^2n)P{^2n ^ Se2n+l) = i4 + i4 = ? ^^^^ 

Similarly it is possible to show 

P{p2n ^ Ce2n) = P{P2n Se2n) = P{p2n ^ Ce2„,-l) = 1/4. (53) 

Thus, after the lapse of time T four levels of the nondegenerate area G will be occupied 
with the identical probabilities 1/4 (Fig. 10.) The motion of the system upwards on energy 
terms will cease upon reaching the level for which the points of the branching in Fig. 10 are 
on the distance at which the condition (34) no longer is valid. The motion of the system 
downwards will be stopped upon reaching the zero level. If the system at the initial moment 
is in the state 2n = N/2, then after N/2 cycles of modulation all levels will be occupied. 
It is easy to calculate a level population for the extremely upper and extremely lower levels. 
Really, the level population for extreme levels is possible to define with the help of a Markov 
chain containing only one possible trajectory in the spectrum of Mathieu characteristics: 

T T 

P{ceN/2, to] Ceiy/2+l, to + CCn, to + N—) = 

T T 
P{seN, to + N-^ ccN-i, to + {N- 1)-)... 

T T 

P{ceN/2+l,to + T 4- seAr/2+l, to + — )^'((se7V/2+l, to + — ^ CeAr/2+1, to), (54) 

where to is an initial time. Here, when discussing the transition probabilities from one state 
to another, we also use a time argument. It is possible to write a similar chain of level 
population for the extremely lower level. As the probabilities of passages, included in the 
right side of (54) by way of factors, are equal to 1/2, then probabilities of an extreme level 
population will be (1/2)^^^. As to the Markovian chain for non-extreme levels, it has a 
cumbersome form and we do not give it here. 



b) Sel f chaotization produced by a big phase "incursion" of the probability amplitude. 
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a) 



FIG. 10: The fragment of the energy terms, participating in passages calculated in the text, a) 
The initial state. The particle is in the state of ce2n- b) The final state. Levels which are affected 
by a change of the field amplitude get populated. 

Different from the area of non-degenerate states G, in the areas of degenerate states G_ 
and G+, the non-diagonal matrix elements of perturbation H'{t) (36) are not zero. 

K- = =< i'L+imm^n+i >~ A/ J ^++i'&iCos2^d^ ^ 0. (55) 

o 

where the wave functions ip± have been defined previously by (30). Here, for the brevity of 
the notation, we omit the upper indices indicating the quantum state. An explicit depen- 
dence of H'{t) on time given by the multiplier cosut is assumed to be slower as compared 
with the period of passages from one degenerate state to another produced by the non- 
diagonal matrix elements H'j^^- Therefore below, perturbations H'(t) will be treated as 
time-independent perturbations able to produce the above-mentioned passages. Therefore 
in the area of degenerate states the system can be found in the time-dependent superposition 



state 
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23j: 



M) = C+(t)V;+ + Q^(t)^2n. (56) 

Probability amplitudes C^{t) are found by means of the following fundamental quantum- 
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mechanical equation expressing the causahty principle: 



-ih^ = {Eon + K+)C^ + 



fin- 

- = H'^_Ct + (i?on + H'__)C-. (57) 

It is easy show that in the case of our problem it should be assumed that H'^j^ = H'__ 
and H'_^_ = H'__^. Let us investigate changes that occurred in the state of the system during 
time AT while the system was in the area G_, (i.e., during the time of movement to the left 
from and, reversal, to the right to /"). It will be assumed that T is part of the period of 
modulation T. For arbitrary initial values the system of equations (57) has a solution 

CM - ^(21±^e.p[^(^ - H'n - - H'n (58, 

where E ^ Eq + H'^^, E ^ Eq + H'__, H'^_ H'. 

After complementing H'{t) (58) with the factor cosz/t, we can take into consideration 
also a slow time dependent change of perturbation (if' H' cosut). 

Let motion begin from the state ■ip2n degenerate area. Then as the initial conditions 

we take 

Q-(0) = 1, C(0) = (59) 
as initial conditions. Substituting (59) into (58), for the amplitudes C±{t) we obtain 

COSCijt, 

C+{t) = -ie('/^)^* sincc^t, uj = ^. (60) 

n 

Now, using (60), for the distribution |V'(^)P (56) we obtain 

|z^(t)r = cos^(^t)|^-(i)r + sin^(f^^)lV^+(i)r- 

-^^n{2^t)[Mt)rjt)-^lj.{t)rM- (61) 

In the expression for |'?/'(t)P the first two terms correspond to the transition probabilities 

> + and — ^ +, respectively, while the third term corresponds to the interference of 

these states. Distribution (61) corresponds to a pure state. 
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Note that (like any other parameter of the problem) the value H' contains a certain small 
error 5H' <^ H', which during the time of one passage 27rh/H' leads to an insignificant 
correction of the phase 27r{SH' / H'). However, during the time T a phase incursion takes 
place and a small error 6H' may lead to uncertainty of phase ~ {6H' /h)AT, which may turn 
out to be of order 2tt. In that case the phase becomes random. Therefore by the moment 
AT the distribution takes the form that can be obtained from (61) by means of averaging 
with respect to the random phase a = {6H' /h)AT. 

Hence, after averaging expression (61), equating the interference term to zero and taking 
into account that 



sm'[{H'/h)t] = cos^[{H'/h)t] = 1/2 

we get 

m' = mw = l{\^+\' + \i^-\% (62) 

where the stroke above denotes the averaging with respect to time. The obtained formula 
(62) is the distribution of a mixed state, which contains probabilities of degenerate states 
\ip±\'^ with the same weights 1/2. The assumption that a large phase is a random value that, 
after averaging, makes the interference term equal to zero is frequently used in analogous 



situations 



23|. 



Thus we conclude that if the system remains in the areas G± of degenerate states for a 
long time, AT ^ 2'n-h/H', AT ^ 2'Kh/5H' during which the system manages to perform 
a great number of passages, then in the case of a passage to the nondegenerate area G 
the choice of continuation of the path becomes ambiguous. In other words, having reached 
the branch point, the system may with the same probability continue the path along two 
possible branches of the Mathieu characteristics. The error H' is evidently connected with 
the error of the modulation amplitude value. It obviously follows that, when passing the 
branch point, the mixed state (62) will transform with a 1/2 probability to the states ce 
and se, as shown in formulas (42) and (43). Analogously, we can prove the validity of all 
subsequent formulas for the passage probabilities (44), (45). 

c) Kinetic stage of evolution 

We have already met with one physical problem that can get reduced to the solution of 
a quantum pendulum (9) - this is a problem of quantum nonlinear resonance. Now we will 
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get to know with other quantum-mechanical problems, that also get reduced to the solution 
of quantum pendulum. 

Note that the parameter cj, which is connected with the modulation depth A/, has (like 
any other parameter) a certain small error 5uj, which during the time of one passage t ~ 
2tx /uo^ leads to an insignificant correction in the phase 2'k{5uj/uj). But during the time 
t ~ AT, there occur a great number of oscillations ( phase incursion takes place) and, in 
the case AT ^ r, a small error 5uj brings to the uncertainty of the phase ~ AT6lj which 
may have order 27r. Then we say that the phase is self-chaotized. 

Let us introduce the density matrix averaged over a small dispersion 6ui: 

.i-wH I . (63) 



where Wj^(t) = |C^(t)P, Fn{t) = \C:i![{t)C~*{t)\. The overline denotes the averaging over 
a small dispersion 5a; 



1 



^(^,t) = ^ / A{x,t)dx (64) 



To solve (64) we can write that 



1 



W+(t) = sin^a;^, W-{t) = cos^o^t, F^t) = -|sin2u;t|. (65) 
After a simple integration of the averaging (64), for the matrix element (65) we obtain 

W^{t) = ^ilTfi26iut) cos2Lut), 

Fnit) = F:{t) = ]^f(25ujt) sin2u;t, (66) 
^, ^ , sin25ujt 

At small values of time t <^ r (r = 2n/6uj), insufficient for self-chaotization {f{25ujt) ~ 1), 
we obtain 

iy+(t < r) = sin^ ut, W:^{t < r ) = cos^ ut, T„(t < r) = ^ sin 2ujt. 

Comparing these values with the initial values (65) of the density matrix elements, we 
see that the averaging procedure (64), as expected, does not affect them. Thus, for small 
times we have 

sin^ ujt ^ sin 2Lut , 
Ptrit«r)=\ ^ ,1- (67) 

sin 2ujt cos cut 
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One can easily verify that matrix (67) satisfies the condition p(t ^ r ) = p{t ^ r), which 
is a necessary and sufficient condition for the density matrix of the pure state. 

For times even smaller than t ^ r ^ r, when passages between degenerate states 
practically fail to occur, by taking the limit ut in (67), we obtain the following relation 
for the density matrix: 



Pr(i = 0) = p+-(t«r)= . (68) 




This relation corresponds to the initial condition (59) when the system is in the eigenstate 
ip2n- Let us now investigate the behavior of the system at times t >t when the system gets 
self-chaotized. 

On relatively large time intervals t > r, in which the self-chaotization of phases takes 
place, for the matrix elements we should use general expressions (66). The substitution of 
these expressions for the matrix elements (66) into the density matrix (63) gives 

/ 1 - f{26ujt)cos2ujt if {26ujt) sin 2ujt \ 

/^n (*) = 2 , , , , • (69) 

\ -if {26ut) sin 2ut 1 + f{26ujt) cos2ut I 

Hence, for times t > T during which the phases get completely chaotized, after passing to 
the limit 6ujt 1 in (69), we obtain 

, [ l-0(e) iO(e) \ , , 

pri»r) = | h 70 

\^ -iO{e) l + 0(£) j 

where 0{e) is an infinitesimal value of order e = 

The state described by the density matrix (70) is a mixture of two quantum states ^/'^ and 
■i/^^ with equal weights. The comparison of the corresponding matrix elements of matrices 
(70) and (69) shows that they differ in the terms that play the role of quickly changing 
fluctuations. When the limit is t ^ r, fluctuations decrease as ~ (see Fig. 11 and 12). 

Thus the system, which at the time moment t = was in the pure state with the wave 
function (68), gets self-chaotized with a lapse of time t ^ T and passes to the mixed 
state (70). 

In other words, at the initial moment the system had a certain definite "order" expressed 
in the form of the density matrix p+~(0) (68). With a lapse of time the system got self- 
chaotized and the fluctuation terms appeared in the density matrix (69). For large times 
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FIG. 11: Time-dependence of the diagonal matrix element W^{t) of the density matrix (63), 
constructed by means of formulas (63), (66) for the parameter values u; = 1/r = 1, C^(0) = 
1, C^(0) = 0. As clearly seen from the Figure, the higher the dispersion value of the parameter 
6uj, the sooner the stationary value W^{t > f 



5uj ■ 



i is achieved. 



t ^ T a, new "order" looking like a macroscopic order is formed, which is defined by matrix 
(70). 

After a half-period the system passes to the area of nondegenerate states G (68). In 
passing through the branch point, there arise nonzero probabilities for passages both to the 
state ce2n and to the state se2n- Both states ip2n ^^^1 ip2n ^i^l contribute to the probability 
that the system will pass to either of the states ce2n and se2n- For the total probability of 
passage to the states ce2n and se2n we obtain respectively 



2lT 



P{ptn (t > ^) ^ Ce2n) = o - / ^tn{'^)ce2n{<^)d(p\ + 

Z IT 



2tt 



1| 1 /" / N / |2 11 11 1 

+^\- hp2ni^)ce2n{^)M =2" 2 + 2' 2 = 2' 



2lT 



1 I 1 / |2 

P{ptn (i > ^) ^ Se2n) = o - / ^tniV^) Se2n{<^)d<^\ + 

Z 71 



2lT 



1| 1 /" / N / |2 11 11 1 

+ ^\- H^2ni^)se2n{^)d^\ = ^ 2 + 2-2 = 2' 



(71) 
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FIG. 12: The vanishing of nondiagonal matrix elements of the density matrix (63) with a lapse of 
time t > f while the system remained in the degenerate area G- . The graph is constructed for the 
parameter values w = l/r = 1, C+(0) = 1, C^(0) = 0., with the aid of formulas (63) and (66). 

Thus, in the nondegenerate area the mixed state is formed, which is defined by the density 
matrix 



where i and k number two levels that correspond to the states ce2n and se2n- 

As follows from (72), at this evolution stage of the system, the populations of two nonde- 
generate levels get equalized. It should be noted that though the direct passage (37) between 
the nondegenerate levels is not prohibited, perturbation (36) essentially influences "indirect" 
passages. Under " indirect" passages we understand a sequence of events consisting a passage 
G through the branch point, a set of passages between degenerate states in the area 

and the reverse passage through the branch point G. The "indirect" passages 

ocurring during the modulation half period T/2 result in the equalization (saturation) of 
two nondegenerate levels. 

As to the nondegenerate area, the role of perturbation H'(t) in it reduces to the displace- 
ment of the system from the left branch point to the right one. 

It is easy to verify that after states (72) pass to the states of the degenerate area G+, we 
obtain the mixed state which involves four states ^fni^) C2n+i(v^)- 

Let us now calculate the probability of four passages from the mixed state C'^O) 




(72) 
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As a result of these passages, in the area we obtain the mixed state described by the 
four-dimensional density matrix 

/l o\ 
10 
10 
\ 1 / 

where the indices of the density matrix (26) show that the respective matrix elements are 
taken with respect to the wave functions ^^(9?) and Ctn+ii^) of degenerate states of the 
area G+. 

It is easy to foresee a further evolution course of the system. At each passage through 
the branch point, the probability that an energy level will get populated is equally divided 
between branched states. We can see the following regularity of the evolution of populations 
for the next time periods. 

After odd half periods, the population of any n — th nondegenerate level is defined as an 
arithmetic mean of its population and the population of the nearest upper level, while after 
even half periods as an arithmetic mean of its population and the nearest lower level. This 
population evolution rule can be represented both in the form of Table 2 and in the form of 
recurrent relations 

P[n, 2k] = P[n + 1, 2k] = ^{P[n, 2A; - 1] + P[n +l,2k- 1]), 



P[n + 1, 2A; + 1] = P[n + 2, 2A; + 1] = -{P[n + 1, 2k] + P[n + 2, 2k]), 



(75) 



where P[n, k] is the population value of the n — th level after time k^, where k is an integer 
number. The creeping of populations among nondegenerate levels is illustrated in Fig. 10. 
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FIG. 13: Results of numerical calculations performed by means of recurrent relations (75). 

This Table is a logical extrapolation of the analytical results obtained in this subsection. It 
shows how the population concentrated initially on one level uq gradually spreads to other 
levels. It is assumed that the extreme upper level rig +4 and the extreme lower level riQ — b are 
forbidden by condition (34) and do not participate in the process. The results of numerical 
calculations by means of formulas (75) are given in Fig. 13 and Fig. 14. Fig. 13 shows the 
distribution of populations of levels after a long time t ^ T when the population creeping 
occurs among levels, the number of which is not restricted by (34). Let us assume that at the 
initial time moment t = 0, only one uq — th level is populated with probability P{nQ) = 1. 
According to the recurrent relations (75), with a lapse of each period T "indirect" passages 
will result in the redistribution of populations among the neighboring levels so that, after a 
lapse of time t = kr ^ T, populations of the extreme levels will decrease according to the 
law which follows from (54) 

If the number of levels defined by condition (34) is finite, then, after a lapse of a long 
time, passages will result in a stationary state in which all levels are populated with 
the same probability equal to 1/A^ (see Fig. 14). Let us summarize the results we have 
obtained above using the notions of statistical physics. After a lapse of time AT, that can 
be called the time of initial chaotization, the investigated closed system (quantum pendulum 
+ variable field) can be considered as a statistical system. 

Formation of statistical distribution of populations of levels P{n) with a lapse of a large 
evolution time t ~ lOOOT of the system. The result shown in this figure corresponds to the 
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TABLE II: Evolution of populations of non-degenerate levels 

case for which the level population creeping is not restricted by condition (34). 

At that, the closed system consists of two subsystems: the classical variable field (36) 
that plays the role of "a thermostat" with an infinitely high temperature and the quantum 
pendulum (7). A weak (indirect) interaction of the subsystems produces passages between 
non-degenerate levels. After a lapse of time t ^ T this interaction ends in a statistical 
equilibrium between the subsystems. As a result, the quantum pendulum subsystem acquires 
the thermostat temperature, which in turn leads to the equalization of level populations. 
The equalization of populations usually called the saturation of passages can be interpreted 
as the acquisition of an infinite temperature by the quantum pendulum subsystem. 

With a lapse of a large time interval t ~ lOOOT the formation of stationary distribution 
of populations among levels takes place. By computer calculations it was found that in the 
stationary state all N levels satisfying condition (34) were populated with equal probability 
1/N. 
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FIG. 14: Results of numerical calculations performed by means of recurrent relations (75). 

G. Analogy between the classical and the quantum consideration. 

The quantum-mechanical investigation of the universal Hamiltonian (mathematical pen- 
dulum), which is reduced to the investigation of the Mathieu-Schrodinger equation, showed 
that on the plane {E,2l) there exist three areas G+, G- and G (see Fig. 8 and 9) differing 
from each other in their quantum properties. Motion in the area of degenerate states G- is a 
quantum analog of rotating motion of the pendulum, while motion in the area of degenerate 
states is an analog of oscillatory motion of the pendulum. The area G lying between G_ 
and G+ can be regarded as a quantum analog of the classical separatrix. The main quantum 
peculiarity of the universal Hamiltonian is the appearance of branching and merging points 
along energy term lines. Branching and merging points define the boundaries between the 
degenerate areas G± and the non-degenerate area G. If the system defined by the universal 
Hamiltonian is perturbed by a slowly changing periodic field, then on the plane {E, 21) the 
influence of this field produces the motion of the system along the Mathieu characteristics. 
If, moreover, the system is in degenerate areas for a sufficiently long time, then the phase 
incursion of wave function phases occurs while the system passes through branching points, 
which leads to the transition from the pure state to the mixed one. As a result of a multiple 
passage through branching points, the populations creep by energy terms (Fig. 10). The 
thus obtained mixed state can be regarded as a quantum analog of the classical stochastic 
layer. The number of levels affected by the irreversible creeping process is defined by the 
amplitude of the slowly changing field. 
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FIG. 15: Analogy between the classical and quantum considerations. Unperturbed motion, a) 
Classical case. Phase plane. Separatrix; b) Quantum case. Specific dependence of the energy 
spectrum on the parameter (Mathieu characteristics). Degenerate G± and non-degenerate G areas 
of the spectrum. 

The classical mathematical pendulum may have two oscillation modes (rotational and 
oscillatory), which on the phase plane are separated by the separatrix (see Fig. 15 a). 

On the plane {E,2l) the quantum pendulum has two areas of degenerate states G- and 
G^. Quantum states from the area G„ possesses translational symmetry in the pendulum 
phase space. These states are analogous to the classical rotational mode. Quantum states 
from the degenerate area possess symmetry with respect to the equilibrium state of the 
pendulum (see axis in Fig. 1) and therefore are analogous to the classical oscillatory state. 
On the plane {E,2l), the area of non-degenerate states G, which lies between the areas 
and G_, contains the line E = 21 corresponding to the classical separatrix (see Fig. 15 
b). If the classical pendulum is subjected to harmonically changing force that perturbs a 
trajectory near to the separatrix, then the perturbed trajectory acquires such a degree of 
complexity that it can be assumed to be a random one. Therefore we say that a stochastic 
motion layer (so-called stochastic layer) is formed in the neighborhood of the separatrix 
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FIG. 16: Analogy between the classical and quantum considerations. Perturbed motion: a) Clas- 
sical case. Stochastic trajectories in the neighborhood of the separatrix form the stochastic layer 
(cross-hatched area); b) Quantum case. The mixed state was formed as a result of population of 
non-degenerate levels situated on both sides of the classical separatrix. 

(see Fig. 16 a). In the case of quantum consideration, the periodic perturbation (36) brings 
about passages between degenerate states. As a result of repeated passages, before passing 
to the area G the system gets self-chaotized, passes from the pure state to the mixed one 
and further evolves irreversibly. While it repeatedly passes through the branch points, the 
redistribution of populations by the energy spectrum takes place. Only the levels whose 
branch points satisfy condition (34), participate in the redistribution of populations (see 
Fig.l6 b). 

H. Investigation of quantum chaos of internal rotation motion in polyatomic 
molecules. 

Let us consider two limiting cases of a low and a high energy barrier. In the limit of 
a low barrier Vo 0, the Mathieu-Schrodinger equation (18) implies the equation for free 
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rotation 

d'^ip 21 , 

+ ^^^^ = °- 

From which for the energy spectrum we obtain Er = (/i^/2/)r^, where r = 0,1,2,... 
are integer numbers. Using the above-given numerical estimates for the molecule of ethane 
C2HQ, we obtain Sr ~ 0.21 ■ 10~^^r^ J, which corresponds to the cyclic frequency of rotation 
^ 2.0 ■ lO^^r'^rad/s. Comparing the expression for energy with the value of the ethane 
molecule barrier we see that only levels with a sufficiently large quantum number (r > 10) 
are located high above the barrier and it is only for such levels that the considered limit is 
valid. 

In another limiting case of a high barrier Vo the rotator is most of the time inside one 
of the potential wells where it performs torsional motions. In that case a can be treated as 
a small angle. After expanding the potential energy in the Schrodinger equation (9) into 
small angles cos 2a ~ 1 — 2a^, we obtain a quantum equation for the oscillator, whose energy 
spectrum has the form Er = {r + l/2)hw, where u = n ^Vo/2I 6.0 ■ lO^^rad/s. For the 
energy spectrum of small torsional oscillations we obtain e,. ^ (r + 1/2) ■ 0.63 ■ 10~^°J. If we 
compare the obtained expression for the spectrun with the corresponding numerical value of 
the barrier, then we can see that only the first two levels Eo ~ 0.32- lO^^'' J, Ei ~ 0.95- lO^^'' J 
and E2 ~ 1.58 ■ lO"^'' J are located in the well but not at a sufficiently large depth that would 
allow us to assume that passages between them correspond to small oscillations. Thus we 
can conclude that for internal rotation of the molecule of ethane C2HQ the approximation of 
small oscillations is not carried out sufficiently well, while the approximation of free rotation 
is carried out for large quantum numbers. 

A real quantitative picture of the internal rotation spectrum can be obtained by means 
of the Mathieu-characteristics when the points of intersection of the line / = lo with the 
Mathieu- characteristics is projected on to the energy axis (see Fig. 17). 

These conclusions are in good agreement with experimental data. In particular, in the 
experiment we observed the infrared absorption by molecules of C2HQ at a frequency ~ 



8.7- 10^^ Hz [24]. For an energy difference between the levels participating in the absorption 
process we have the estimate AEexper ~ 0.54 ■ 10~^°J. Comparing the experimental result 
with the energy difference between two neighboring levels in the case of approximation of 
small oscillations, we obtain 

Aeo = 2ni'oh, 
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Inserting the parameters values for molecules of C2HQ, we obtain the estimate Aeo = 



0.63- lO-^V. 

After comparing the obtained estimates with the energy difference between the states 
described by the wave functions ce^lipjo), se^lipjo) and applying formula (20), we obtain 
Aeicesi^, Q < > se^i^, Q) ~ 0.54 ■ lO-^o J. 

Comparative analysis of the obtained estimates gives us the grounds to conclude that 
the energy levels corresponding to the states ce^i^ip, lo) , se^^ip, lo) participate in the infrared 
absorption revealed in the experiment. 

As we see from Fig. 17 energy levels involved in the process are distributed in a random 
manner. In some sense, this resembles RMT. However as number of levels is small, and 
energy spectrum is well defined. 
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Let us assume that the considered quantum system is subjected to a radio frequency (RF) 
monochromatic pumping whose frequency Q satisfies the condition Q <^ Vo/h. This causes 
a slow modulation of quick electron motions in a molecule. The formation of an energy 
barrier is a result of the averaging over quick electron motions and thus it is obvious that 
due to the pumping effect the barrier value is time-dependent, 

Vo^Vo + AVcosnt. (76) 

The depth of modulation AV depends on a pumping power. By replacing (76) we obtain 
the time-dependent Hamiltonian 



^oiv) = + lo COS 2y?, (77) 



dip 

H'{(f,t) = A/ cos 2(p cos fit. 



Al = -^AV. (78) 

Simple calculations show that the matrix elements of perturbation H'{ip,t) with respect 
to the wave functions of the nondegenerate area G are equal to zero 

{cen\H'{ip,t)\sen) ^ Al J ceni^p) cos2ipsen{ip)d(p = 0, C^^) 

o 

where n is any integer number. Therefore perturbation (78) cannot bring about passages 
between nondegenerate levels. 

The interaction H'{(f,t), not producing passages between levels, should be inserted in 
the unperturbed part of the Hamiltonian. The Hamiltonian obtained in this manner can be 
considered as slowly depending on time. 

Thus, in the nondegenerate area the Hamiltonian can be written in the form 

H = —^ + l{t) cos 2<f, (80) 

Because of the modulation of the parameter l{t) the system passes from one area to 
another, getting over the branch points. 

As different from the nondegenerate state area G, in the areas of degenerate states 
and the nondiagonal matrix elements of perturbation H'{t) (78) are not equal to zero. 
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For example, if we take the matrix elements with respect to the wave functions ipfn+i 5 then 
for the left degenerate area G- it can be shown that 

K- = H'.+ = (^2Vll^'(¥'>^)l^2n+l) ~ A^/ i^tn+li^2:+lCOs2^d^ ^ 0. (81) 

o 

Note that the value H'^^ has order equal to the pumping modulation (78) depth Al. 

Analogously to (81), we can write an expression for even 2n states as well. 

An explicit dependence of H'{(f,t) on time given by the factor cosQt is assumed to be 
slow as compared with the period of passages between degenerate states that are produced 
by the nondiagonal matrix elements -f^+_- Therefore below the perturbation H'_^__{ip,t) will 
be assumed to be the time-independent perturbation that can bring about passages between 
degenerate states. 

In a degenerate area the system may be in the time-dependent superpositional state 

V^2n(t) = C^mtn + (i)fe- (82) 

The probability amplitudes C^{t) are defined by means of the fundamental quantum- 
mechanical equation, expressing the casuality principle. We write such equations for a 
pair of doubly degenerate states: 

- ih^ = H'^_C^ + (i?,„ + i/:„)C-. (83) 

where the matrix elements are taken with respect to degenerate wave functions (see (81)) 
and Eon is the energy of the n — th degenerate level near a branch point. 

Let us assume that at the initial moment of time the system was in the degenerate state 
il)2n- Then as initial conditions we should take 

C„"(0) = 1, C+(0) = 0. (84) 

Having substituted (82) into (81), for the amphtudes we obtain 

C^{t) = iexp{—Et) smut, 

C-{t) = expi^Et) coscjt, (85) 
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E = Eon + H', 



where u = ^ = is the frequency of passages between degenerate states, r is the passage 
time. 

Note that the parameter u has (hke any other parameter) a certain small error Su, 
which during the time of one passage t ~ 211/00, leads to an insignificant correction in the 
phase 2'k{5uj/uj). However, if during the time t ~ AT, when the system is the degenerate 
area (AT < T, T = 27r/f2) there occurs a great number of passages (AT ^ r), then for 
AT5uj ~ 27r, a small error 5uj leads to the phase uncertainty. Then we say that the phase is 
self-chaotized. The self-chaotization formed in this manner can be regarded as the embryo 
of a quantum chaos which, as we will see in the sequel, further spreads to other states. 

After a half-period, the system passes to the area of nondegenerate states G. In passing 
through the branch point, there arise nonzero probabilities for passages both to the state 
ce2n(v^) and to the state se2n(v^)- Thus, in the nondegenerate area the mixed state is formed, 
which is defined by the density matrix 



where i and k number two levels that correspond to the states ce2n(v^) and se2n(v^)- 

As follows from (86), at this evolution stage of the system, the populations of two non- 
degenerate levels get equalized. It should be noted that though the direct passage (79) 
between the nondegenerate levels is prohibited, perturbation (78) essentially influences "in- 
direct" passages. Under "indirect" passages we understand a sequence of events consisting 
a passa ge G ^ through the branch point, a set of passages between degenerate states 
in the area (?_, and the reverse passage through the branch point G^ G. The "indi- 
rect" passages occurring during the modulation half period T/2 result in the equalization 
(saturation) of two nondegenerate levels. 

Thus "indirect" passages are directly connected with a quantum chaos. Hence, by fixing 
"indirect" passages we thereby fix the presence of a quantum chaos. 

Let us assume that the investigated molecule is a component of a substance in a gaseous 
or liquid state. Then the molecular thermal motion, which tries to establish an equilibrium 
distribution of populations according to Boltzman's law, will be a "competing" process for 




(86) 
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FIG. 18: A thermodynamic scheme of the process. Subsystem 1 is a usual thermostat with tem- 
perature Tq , subsystem 2 is a thermostat having an infinite temperature and consisting of the 
interaction H'{ip,t) (78); the subsystem 3 is the quantum system corresponding to internal rota- 
tions of molecules and being able to receive energy from subsystem 2 and to transfer it to subsystem 
1. 

the quantum chaos described above. Using thermodynamic terminology, we can say that the 
considered quantum system is located between two thermostats. One of them with medium 
temperature Tq tries to retain thermal equilibrium in the system, while the other, having an 
infinite temperature, tries to equalize the populations. 

An equation describing the change of populations according to the scheme shown in Fig. 18 
has the form 

^ — =Sf^. (-) 

where 2W is a probability of "indirect" passages and Ti is the time of thermal chao- 



tization. Bloembergen, Parcell and Pound used equation 



saturation of nuclear magnetic resonance in solid bodies 25||. For a stationary distribution 



87) to describe the process of 
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of populations from (87) we obtain 

n. = nf)^, (88) 

where s = 2WTi is called the saturation parameter. For s ^ 1 "indirect" passages have a 
stronger effect on the system than thermal processes. 

Difference from the magnetic resonance consists in the expression of the transition prob- 
ability W. Whereas, in the case of magnetic resonance, transition probability 2W ~ Hf/A 
is proportional to the pumping intensity Hi and inversely proportional to the line width A, 
in the case of non-direct transitions, transition probability 2W ~ 1/T does not depend on 
the oscillation amplitude and is inversely proportional to the oscillation period T (i.e. 
to the time of motion between branch points). But, it is worth to keep in mind, that the 
amplitude of oscillation must be large enough for the system to achieve the branch points. 

Thus, along with the conditions that the perturbation is adiabatic Q <^ Vo/h, phase 
incursion AT ^ r and self-chaotization during the process of multiple transitions between 
the degenerated states takes place AT6u! ~ 27r, 6u! <^ u, AT < T. the condition s ^ 1 is 
a necessary condition for the formation of a quantum chaos. In the opposite limiting case 
s <^ 1, nj ~ nfK the quantum chaos will be completely suppressed by thermal motion. 
Let us proceed to discussing the quantum chaos possible experimental observation. The 
most suitable material for this purpose in our opinion is ethane C2HQ. But, by use of 
easy estimations one can prove, that for the ethane in the gaseous state, conditions of the 
quantum chaos observation hardly can be achieved. Really the requirement that modulation 
must be slow Q ^ Vq/H ^ lO^^sec"^, leads to the necessity to take the modulation frequency 
from the transition radio frequency range ^ 0.3 ■ lO^^sec"^. 

In the case of gases the thermal chaotization time is estimated by the formula Ti ^ ^ 
, where d is the molecule size, u = Vi^ = A/3/2(2A;To/m)^/^ is a mean motion ve- 
locity of molecules. After substituting the numerical values for C2HQ, we obtain Ti ^ 
0.3 ■ lO^^VV^o, ^0 ~ 200°K > 184" K, (where 184°i^ is the boihng point of the ethane at 
normal pressure). Then for the time of thermal chaotization and for the saturation param- 
eter we get respectively: Ti ^ 0.3 ■ 10~^^sec, S ~ (Ti/T) ~ 1. Hence, thermal chaotization 
and quantum chaos are equally manifested in the system. 

In the case of a liquid under Ti we should understand the mean time of the settled life of 
a molecule, which is about 10~®sec. Relatively large times of relaxation in liquids ensure the 
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fulfillment of the saturation condition (s ^ Ti/T ^ 10^). In such a way, for ethane C2HQ in 
the liquid state, quantum chaos caused by the non-direct transitions has stronger influence 
upon the system then common thermal chaotization. At the same time, the state with the 
equally populated levels will be formed in the system. 

Suppose, the system (ethane in the liquid state) is subjected to the action of a infrared 
radiation field and infrared absorption at the frequency ~ 289sec~^ is observed 2J]. This 
absorption corresponds to the transitions between the states se-^ <-> ces. Furthermore, let 
us suppose that along with the infrared, the radio frequency pumping which can cause the 
non-direct transitions is applied. This in its turn leads to the discontinuance of the infrared 
absorption. This phenomenon can be considered as an observation of quantum chaos. 

Now we consider the example of the reaction of isomerization LiCH CHLi, that was 
discussed in subsection E. As was mentioned above the most simplified model, corresponding 
to this process, can be presented with the aid of Mathieu-Schrodinger equation (77). It is 
known, that energy barrier of transition of the process of vibration of the atom Li in the 
process of isomerization \^ ^ 7 - 10^^^ J. The moment of inertia is Jlj = 1.1 ■ lO'^^mcg"^ That 
is why, in the case = y^^VVLi ~ 19 number of vibration energy levels are packed in the 
potential well. Now it is not hard to determine energy spectrum with the aid of Mathieu 
characteristics. 

For this process the presence of chaotically distributed levels in infrared region was es- 
tablished with the aid of numerical experiments. 

Now let us deal with the problem of nonlinear quantum resonance from which we have 
begun review. To begin with let us estimate resonance value of Jq, close to which variation of 
values of AJ may take place. By definition resonance value of action Iq is determined from 
equation uj — ujq + uJn{Io) = 0. By taking typical values of parameters, know from optics: 
uo ~ 10^^ sec~\ u -uo 10^^ sec-\ 7 ~ 4 ■ 10^^ J /kg ■ we obtain: Iq = ~ 
0.5 ■ 10^^^ J ■ sec. This value of action I/H c:^ 0.5 ■ 10^ times as much as Plank constant. So, 
resonance takes place on very high levels (in quasi-classical region) of atoms potential well. 
As we have done in previous cases, it is possible to find variation spectrum of action AI 
with the aid of Mathieu characteristics. For this, we must know the value of dimensionless 
parameter Iq = 2V (Iq) / u' h"^ . It is easy to find the value tu' ~ 2 ■ 10^^ J~^sec~^ and for strong 
light fields /o — lO^V/m. As a result we obtain Iq ~ 10^°. This means in its turn, that 
there are a large number of levels in potential well. Because of the large value of Iq, it is 
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FIG. 19: Energy spectrum of the process of isomerization LiCH CHLi , obtained as spectrum 
of eigenvalues of Mathieu-Schrodinger equation, when Iq = -^Vq ~ 350. Speed of energy levels 
close to E = 21q 700 can easily be observed in figure. 

not possible to determine spectrum by method, used previously, with the aid of Mathieu 
characteristics. However, by extrapolating it is possible to expect, that at approaching the 
top of well, alongwith the common attraction of levels, their chaotization will take place. In 
our opinion, chaotically distributed levels of great density are the basis for random matrix 
assumption. 

II. THE PECULIARITIES OF ENERGY SPECTRUM OF QUANTUM CHAOTIC 
SYSTEMS. RANDOM MATRIX THEORY. 

In the first part of this work we have considered the possibility of formation of mixed 
state in quantum chaotic system. The peculiarities of spectral characteristics of the univer- 
sal Hamiltonian (7) of concrete physical system corresponding to the Mathieu-Schrodinger 
equation (9) was a cause of this. 

In the present part we shall try to get the same results proceeding from more general 
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consideration. In particular, we shall use methods of random matrixes th eory [?, |8|, 26|. 



Study of quantum reversibility and motion stability is of great interest [27|] . This interest 
is due to not only the fundamental problem of irreversibility in quantum dynamics, but also 
to practical application. In particular, it reveals itself in relation to the field of quantum 
computing 
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A quantity of central importance which has been on the focus of many studies 29|, |30|, 



32| is the so-called fidelity f{t), which measures the accuracy to which a quantum state 



can be recovered by inverting, at time t, the dynamics with a perturbed Hamiltonian: 

fit) = I < 4!\e'^'e-'^°'\^ > |2, (89) 

where ip is the initial state which evolves in time t with the Hamiltonian Hq while H = 
Hq + V is perturbed Hamiltonian. The analysis of this quantity has shown that under 
some restrictions, the series taken from f{t) is exponential with a rate given by the classical 



Lyapunov exponent [26|]. But here a question appears. The point is that the origin of the 
dynamic stochasticity, which is a reason of irreversibility in classical case, is directly related 
to the nonlinearity of equations of motion. For classical chaotic system this nonlinearity 
leads to the repulsion of phase trajectories at a sufficiently quick rate l2|. 

In case of quantum consideration, the dynamics of a system is described by a wave 
function that obeys a linear equation and the notion of a trajectory is not used at all. 
Hence, at first sight it seems problematic to find out the quantum properties of systems, 
whose classical consideration reveals their dynamic stochasticity. 

In this part of the paper, by using of method of random matrix theory (RMT), we try 
to show that quantum chaotic dynamics is characterized by the transition from a pure 
quantum-mechanical state into the mixed one. 

With this purpose we shall consider a case when the system's Hamiltonian may be pre- 
sented in the form: 

H{t) = Ho + V{t), (90) 

where Hq is chaotic Hamiltonian with irregular spectrum, V^(t + To) = V{t) periodic in time 
perturbation. We shall try to show that in this case irreversibility in the system appears 
as a result of loss of information about the phase factor of the wave function. At the same 
time, unlike of the second part of the paper where we have proved this fact proceeding from 
the spectral peculiarities of the concrete system, now we shall prove this in more general 
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manner. For this reason we shall make use of the fact that according to RMT the eigenvalues 
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34|. 



of chaotic Hamiltonian can be considered as a set of random number 

After taking (90) into account, the solution of the time-dependent Schrodinger equation 

= mm) > (91) 

can be written formally with the help of a time- dependent exponential 

U{t) = exp[-i [ dt'H{t')], (92) 
where the positive time ordering requires: 







, A(t)B(t') if t>t' 
[A{t)B{t%={ 93 
B{t')A{t) if t<t' 

In our case H{t + T^k) = H{t), k = 1,2, .. . the evolution operator referring to one 
period Tq, the so-called Flouqet operator U{Tq) = F is worthy of consideration, since it 
yields the stroboscopic view of the dynamics 

mkTo) >= {py^im > . (94) 

The Flouqet operator being unitary has unimodular eigenvalues. Suppose we can find eigen- 
vectors > of the Flouqet operator 

< V^nlv^m >= Snm- (95) 

Then, with the eigenvalue problem solved, the stroboscopic dynamics may be written ex- 
plicitly 

li^ikTo) >= Y.e-'"^- < ^n\m\Vn > . (96) 

n 

As it was mentioned above, our aim is to prove that one of the signs of the emergence 
of quantum chaos is a formation of the mixed state. Being initially in a pure quantum- 
mechanical state, described by the wave function \ip >, the system during the evolution 
makes an irreversible transition to the mixed state. 

The information about whether the system is in the mixed state or in the pure one, may 
be obtained from the form of the density matrix 3^. Using (96) as a density matrix of the 
system we get the following expression: 

Pnm{t = kTo) = A^^e"'^'^--'^-\ (97) 
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^nm =< V'ml^(O) >< V(0)|</^n > • (98) 

Exponential phase factors of the non-diagonal matrix elements Pnm{t) express the prin- 
ciple of quantum coherence {sgI, and correspond to the complete quantum-mechanical de- 
scription of the system in pure quantum-mechanical state. While they are not equal to zero, 
the system is in the pure state. So, to prove the formation of the mixed state one has to 
show zeroing of non-diagonal elements of density matrix. 

According to the main hypothesis of the random matrix theory [3], the phase /(n, m) = 
'^m — in the exponential factors of the non-diagonal matrix elements in (97) is a random 
quantity. So, it is clear that the values of matrix elements of density matrix of the chaotic 
quantum-mechanical system Pnm{k) are random values too. Taking a statistical average of 
expression (97), we have 

< Pnmik) > = < A^mC''^'^--''-^ > . (99) 

Then in case of random phase f{n,m) has a normal dispersion, we get [s^]: 

Pnm(A;)~e-^. (100) 



This phenomenon is connected with the "phase incursion". Uncertainty of phase in (97) 
is accumulated little by little with time. Finally, when k > V^/cr, the uncertainty of phase 
is of order 27r , and the phase is completely chaotized. As a result the system passes into 
mixed state. 

According to the theory presented above, calculations were made for concrete physical 
„ I„ pa..c..a. in wo.. Q cKaotic syste. of connected osciHato. was studied, and 



in work [38| Kepler's asymmetry problem. The numerical results obtained in these works 
prove the formula (100) to be correct. 

After formation of the mixed state the quantum-mechanical consideration loses its mean- 
ing and there is a need to use a kinetic description. 

For derivation of master equation let us split the density matrix operator p on a slow p^i 
and a fast varying pj^^ operators: 

P = PR + PNR- (101) 

Relevant part pn in the basis \(f)n) contains only diagonal elements, whereas nonrelevant 
part pnr contains only nondiagonal elements. These elements, as it was shown in previous 
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section, contain fast oscillating exponents and when taking average over the ensemble, the 
zeroing of them takes place. Elimination of the diagonal part from the density matrix is 



linear operation, which satisfies the property of projection operator = D 37j 



PB = Dp, pMR = {l-D)p. (102) 

Let us note that this reflection is nonreversible. Due to the zeroing of nondiagonal part 
of the density matrix a part of information is lost. 

Inasmuch as relevant statistical operator pR{t) is different from the total operator p(t), 
generally speaking it does not suit Liouville-Fon Neumann equation 

% + ^Lp = 0, (103) 

where L is Liouville operator. After acting on the equation (103) with D operator, we get 

dpR 



35| 



dt 

dpNR , 



+ iDL{pR + pnr) = 0, (104) 
i{l-b)L{pR + pNR) = Q. (105) 



dt 

For the purpose to obtain closed equation for pR we exclude from the equation (104) p^R 
As a result we get 



dpRjt) 
dt 



+ iDLpR{t)+ [ K{t-ti)pR{ti)dti = -iDLexp[-'i{t-to){l-D)L]pr,R{to), (106) 

J to 



Kit - ti) = DLexp[-i{t - ti){l - D)L]{1 - D)L. (107) 

This equation is valid for t = > r^. Evidently pR{t) is expressed by way of values of 

pR{ti) taken for the time interval to < ti < t, and additionally through the value of pNR{to)- 
If in the initial moment of time t = to the system is in a pure quantum-mechanical state, 
then pnr{Q) ^ 0. 

In this case solving the equation (107) is problematical. However let us recollect that for 
to > V^/cr system is already in a mixed state. Therefore equation (107) takes more simple 
form (p^RiO) = 0) 

+ tDLpRit) + [ Kit- t^)pRit^)dt^ = 0. (108) 



to 



For solving equation (108) we shall use a method of superoperators [39|. But before per- 



forming this, we should note, that we have obtained a closed equation for relevant part of 
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statistical operator. We were able to come to this only because when t > to = \/2/a the 
system is in the mixed state and all nondiagonal matrix elements of the density matrix are 
equal to zero. 

Further when studying the evolution of the system we shall consider as the origin of time 
the moment of the formation of the mixed state in the system. This corresponds to a formal 
transition to limit to — oo. Further for simplification of (108) we shall use Abel theorem 



hm 1 r fm = f{0)- hm r e<'^^f'{t')dt\ 



(109) 



Taking (109) into account, equation (108) will have the following form: 

+ iDLpnit) = lim T e'^''-''^K{t - t')pn{t')dt' . (110) 

dt e^+co 

According to the method of superoperators [4^, the correspondence of one operator to 
another may be considered as representation. The operator in this case will be represented 
by a matrix element with two indices, while the linear product of operators is a matrix with 



four indices, i.e. a superoperator 40|]. A concrete example is the projection of D operator 
on diagonal elements. We should notice, that in our case, the projection operator is some 
definite physical procedure of averaging matrix elements of the density matrix over Gaussian 
chaotic ensemble. 

From the relation Dpnrn{t) = Pnrn{t)5mn we come to the following representation of D 

superoperator Dnmn'm' = ^nn'^mm'^nn SO that 



^ ^ Dmnm'n' Pm'n' Pnm^nm- (-^-^-^) 



Taking (111) into account, (110) is 

+ zD{Lpn)nn = - f_ dt'[K{-t')pn{t + t')Ue''' . (112) 

In expression (112) and further for short we shall omit R index for diagonal matrix elements 
of the operator p. Considering the relation 

[L, PR]nn = LnnaaPaa = (113) 

a 

and representing Liouville operator as L = Lo + L' in compliance with (113) we get 

- f dtie'^'-'''>J2Knnmmit-tl)pmnitl), (114) 
J — oo „ 



dpnnit) 

dt 
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and 

Knnmmif) = [L 6 ^ (1 — D)L ]nnmm- (H^) 

From (115), (114) and from representation of Liouville operator in the form L = + L' , 
one can see, that kernel Knnmmif) is at least of second order by Axq. Next it is easy to 
check correctness of the expression 

L'ahmm = ^^(Km^fem - VmbSam) = Vab - Kb = 0, (116) 
m m 

for Liouville superoperator 

The relation (116) in its turn leads to the rule of sums J2m Kmmnn = 0. Taking the symme- 



tries Labcd = Lcdab, Dabcd = D^dab mto accouut, from (76) we get [37| 

^^^^ = - j rftie^(*-*i) ^[Knnmm{t - tl)p,„„.(tl) " Kmmnnit - tl)pn„(tl)]. (118) 

Next we shall make the following approximations. With accuracy up to the value of (Axq)^ 
order in the exponent in expression (115), we shall replace the complete Liouville operator 
L = Lq + L' with Lq. With the same precision we may set pnn(t — t') = Pnn(t). As a result 
from (118) we get 

dpnnit) ^ J2[W^^p^^^^t^ _ W^^pnnit)], (119) 



dt 

va^n 



where 



W^^ = - r rftV*'[LV*(^^^)^«(l - D)L'Umm. (120) 

J —oo 

Since 

{DLo)abcd = Sabi^aSacSbd - ^b^bdSac] = 0, (121) 

in expression (120) the operator D in the argument of exponential function may 
be omitted. Taking into account the time dependence of the operator x{t) = 



XoAxofit), f{t) Yli e*"^^*, VL = 2'k/To, for nondiagonal matrix elements we get [37|: 

= [WnmPmmit) ' WmnPnn{t)], (122) 

where 

oo 
1/=— oo 
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is the transition amplitude between the eigenstates of the Hamiltonian Ho-, Enm = '^n — 
'^m, Vnm IS the matrix element of the operator Vo = AxoQiQl the basic of eigenfunctions 
of the Hamiltonian Ho, Vnm = {i^n\Vo\i^m) ■ 

Equation (122) describes a nonreversible evolution of the system from nonstationary state 
to the stationary state defined by the principle of detail equilibrium. To prove irreversibility 



of the process let us consider time dependence of nonequilibrium entropy 



39] 



S{t) = -KBj2Pnn{t)ln{pnn{t)), (124) 

n 

where Kb in the Boltzmann constant. Taking into account 'YlinPnn{t) = 1, from (124) we 
get 

dS{t) 



dt 



-Kb Wnm[pmm{t) " Pnn{t)]ln{pnn{t)) - 



-KbJ2^^^ = lKBj2J2^nm[Pnn{t)-pmm{t)][ln{pnn{t))-ln{pmm{m. (125) 

n n m 

Due to the property of logarithmic function 

{Pnn(t) - Pmm(t)){ln{pnn(t)) - ln{pmm{t))) > 0, (126) 

we see that 

^ > (127) 

dt - ^ ' 

This testifies the growth of entropy during the evolution process. 

More exact estimation of entropy growth may be obtained from the principle of detail 
equilibrium WnmPmm = WmnPnn- Taking into account that in our case Wnm = Wmn, for the 
entropy growth we get 



AS = S{t > V2/a) - S{t = 0) = KbIuN, (128) 

where is the number of levels included into the process. 

Thus, we can conclude, that the complex structure of energy spectrum of quantum chaotic 
system leads to the fact that after averaging over small spread in values of parameters there 
occurs formation of a mixed sate in the system. 
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III. CONCLUSION 



The traditional notion of an area, where the laws of statistical physics are effective, 
consists of the assumption that the number of interacting particles is sufficiently large. 
However, a lot of examples of systems with a small number of degrees of freedom, where 
chaotic motions occur, had become known by the end of the last century. A new stage 
in the development of notions about chaos and its origin appeared in the last two decades 
of the last century. It turned out that the classical Hamiltonian system may experience 
a special kind of instability. Because of this instability various dynamic characteristics 
of the system randomly change with time. Such a property of the system that performs 
random motion is called dynamic stochasticity. It is well known that the appearance of 
non-reversibility in classical chaotic systems is connected with the local instability of phase 
trajectories relatively to a small change of initial conditions and parameters of the system. 
Classical chaotic systems reveal an exponential sensitivity to these changes. This leads to 
an exponential growth of the initial error with time, and as the result after the statistical 
averaging over this error, the dynamics of the system becomes non-reversible. In spite of this, 
the question about the origin of non-reversibility in quantum case remains actual. The point 
is that the classical notion of instability of phase trajectories loses its sense during quantum 
consideration. Therefore more detailed analyzes of the possible mechanisms of irreversibility 
in chaotic quantum-mechanical systems is needed. More over even the fact of possibility of 
emergence of irreversibility in chaotic quantum-mechanical systems is not evident. Quantum 
dynamics is unitary and therefore reversible. But on the other hand, without assumption 
of small initial dispersion of the system's parameters, classical chaotic dynamics also will 
be reversible. Note that while studying classical chaos, one usually examines the stability 
of the system relatively to a small change of initial conditions and system parameters. A 
small initial error of these parameters always exists and remains unavoidable (one may 
measure the parameters of the system at a very high precision, but even in this case there 
is still a small error, the removal of which, i.e. the measurement at absolute precision, is 
impossible). So not the existence of the unavoidable error is fundamental, but what kind 
of influence it brings over the system dynamics. It is well known that in case of regular 
systems such an influence is negligible, but if the system is chaotic, the effect of it increases 
exponentially. In this review paper we touched upon results obtained in our previous works. 
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We studied paradigmatic model for quantum chaos, time dependent universal Hamiltonian 
(Hamiltonian of mathematical pendulum). Main result obtained in those works is that reason 
of irreversibility in time dependent quantum dynamics of the chaotic quantum-mechanical 
systems, is a very specific spectral characteristic of the system. For universal Hamiltonian 
and corresponding to it Mathieu-Schrodinger equation above mentioned means very specific 
dependence of the energy eigenvalues on the values of potential barrier (see. Fig.l). In 
addition, transition form pure quantum-mechanical state to the mixed one, was found to be 
quantum analog of the classical chaotic layer, characterized classical chaotic motion near the 
separatrix. Results obtained are important not only for better understanding of physics of 
quantum chaos, but for practical implementation to organic molecules, quantum computing 
and so on. 
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